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Abstract. Using the Lorenz equations as an example we 
show that one chaotic system can be controlled by synchro¬ 
nizing its behavior with the chaotic behavior of another 
system. We particularly discuss the implications of this 
phenomenon in geophysical systems. 


The difficulty of carrying out long-term predictions of 
atmospheric dynamics and the evolution of climate is a 
problem of obvious concern. Nowadays there is increasing 
awareness that deterministic chaos might provide a possible 
paradigm for the complexity of atmospheric and climatic 
dynamics. Periodicity is not the first apparent characteristic 
of the behavior of many geophysical fluid dynamic systems, 
but atmospheric and oceanic flows often exhibit substantial 
coherent features, localized in either or both of space and 
time, which occur sporadically and unpredictably but with 
a certain statistical regularity which can be important in 
extended-range atmospheric prediction. Such features are 
exemplified by blocking patterns in the mid-latitude atmo¬ 
sphere, or by persistent anomalies of the ocean-atmosphere 
system (of which El Nino is the spectacular example 
[McCreary and Andeson, 1991, Brindley et al., 1992]); their 
presence coincides with temporary and localized improve¬ 
ment in potential predictability. 

In this letter we propose a mechanism for reduction in 
chaos which could affect atmospheric potential predictability 
based on the continuous chaos control scheme [Pyragas, 
1992, 1993, Qu et al., 1993]. 

We consider two chaotic systems, which we call A and B 
respectively, 


x = f(x) 

y - g 00 


(l) 


where x,yeR n , and we use the controlling strategy which is 
schematically illustrated in Figure 1; the two systems are 
coupled through the operators X, /*, which we take to have 
a very simple linear form. We assume that some or all state 
variables of both systems A and B can be measured, so that 
we can measure signal x(t) from the system A and signal y(t) 


from B, and that the systems are coupled in such a way that 
the differences D 12 (0 between the signals x(t) and y(t ) are 
used as control signals 

F,(f) = k[x(t) -y(t)] = XD x (f) 

( 2 ) 

F 2 {t) = (i[y(r) - x(t) ] = \iD 2 (t) 

introduced respectively into each of the chaotic systems A 
and B as a negative feedback. We take X, /x> 0 to be 
experimentally adjustable weightings of the perturbation. 

Using the coupling schematically shown in Figure 1 we 
have shown that one chaotic system coupled with the other 
one can significantly change the behavior of one of them 
(unidirectional coupling, i.e., X or /*=0) or of both systems 
(mutual coupling, i.e., X, 0). This property allows us to 
describe the above procedure as the controlling chaos by 
chaos method. 

Propositions 1 and 2 presented in the Appendix give 
rigorous conditions under which chaotic attractors of systems 
A and B are equivalent (Proposition 1), or the evolution of 
one of them is forced to take place on the attractor of the 
other one (Proposition 2). Detailed investigation of the 
question of equivalence of chaotic attractors is given 
elsewhere [Kocarev and Kapitaniak, 1994]. Here we 
describe some applications of controlling chaos by chaos in 
geophysical systems, pointing out that, even when the 
conditions of Propositions 1 and 2 are not fulfilled the 
introduction of coupling can still have practical importance. 

In our numerical examples we first consider two Lorenz 
models [Lorenz, 1963, 1965] mutually coupled in the 
following way 

X x = -aX x + oY l + X(X 2 - X x ) 

#i = - X l - Z 1+ r j X, - Y, + k(Y 2 - YJ 
Z t = Xj Yj - bZ t + X(Z 2 - Zj) 

(3) 

X 2 = - aX 2 + oY 2 + |i (X t - Xj) 

Y 2 = -X 2 - Z 2 + r 2 Xj + |x(Y x - Yj) 
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Z 2 = X 2 - Y 2 - bZ 2 + |i (Zj - Z 2 ) 

where a, r 12 and b are constants. The Lorenz model has 
often been proposed as a paradigm for the "chaotic" extra- 
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Figure 1. Scheme of controlling procedure. 

tropical atmospheric circulation [Palmer, 1993]. The 
variables X, Y and Z then represent in some broad sense 
Rossby wave components of the extratropical general 
circulation. Coupling between two Lorenz models introduced 
in eq. (3) might then be interpreted as mutual interdepen¬ 
dence of extratropical circulations in two regions character¬ 
ized by different r parameter value, say an intensive storm 
track and a relatively stable anticyclonic region. The concept 
of teleconnections of this kind, achieved through the mecha¬ 
nism of quasi-linear Rossby trains, has both theoretical and 
observational support [Madden and Julian, 1971, Weick- 
mann, 1991]. 

Numerical computations have been carried out using 
software INSITE [Parker and Chua, 1989]. In Figure 2 (a-b) 
we show the chaotic attractors of single Lorenz models 
(X,/x=0) for <t= 10.0, ^ = 197.4, 6=8/3 (Figure 2(a)) and 
r 2 =211.0 (Figure 2(b)). These attractors are characterized 
by the following spectra of Lyapunov exponents Xj=1.87, 
\2=0, X 3 = -15.54 (Figure 2(a)) and X^O.78, X2=0,X 3 = 
-14.44. In Figure 2(c) we show the behavior of both above 
mentioned Lorenz systems coupled with X= 100 and fx= 1. 
Although this attractor is still chaotic (^=0.79, X^O, X 3 = 
—14.34), trajectory behavior on it is more predictable as its 
Lyapunov dimension, d L =2.053, is smaller than the dimen¬ 
sion the of original attractor (d L =2.121). This dimension in¬ 
crease is produced by a significant decrease of positive 
Lyapunov exponent (Xj=0.79 in comparison with X^l.87 
of the original attractor). 

In a second example we consider the coupling of a Lorenz 
system with a linear oscillator 

W = -Q (t)V - /t(w - W m ) 

(4) 

v = n(?)(w - w) - kv 

which in geophysical context represents the tropical atmo¬ 
sphere [Palmer, 1993, Madden and Julian, 1971, Weick- 
mann, 1991]. Here, Q is taken to be the frequency of some 
dominant internal mode of large-scale variability of the 
tropics, e.g., the Madden-Julian oscillations. In our compu¬ 
tations we consider 0 to be a time dependent random 
variable with uniform distribution in the interval [1.3, 1.7] 
in nondimensional time. V and W represent two phase- 
quadrature components of the tropical oscillations. For 
example W can be considered as representing the Walker 
circulation [Palmer, 1993]. 

Although eq. (4) is stochastic, any particular phase space 
trajectory has the properties of a chaotic trajectory, so we 



Figure 2. Evolution of Lorenz system; (a) unmodifird first 
system: a= 10.0, r= 197.4, 6=8/3; (b) unmodified second 
system: <r=10.0, r=212.0, 6=8/3; (c) attractor of the first 
system controlled by the second system: X= 100.0, ^=1.0. 

can apply system (4) to control the chaotic behavior of a 
Lorenz model, in the same way as before. 

Application of our chaos control method requires consid¬ 
eration of the following coupled equations: 

X = - aX + aY 

Y = -XZ + rX - Y + \ {W - Y) 

Z = XY - bZ (5) 

W = -QV - k{W - W*) 

V = Q(f) - (W - W*) - kV 

Examples of numerical calculations for <7=10.0, r= 197.5, 
6=8/3 and W* =0 are shown in Figure 3. Previously in 
Figure 2(a) we showed the original attractor of a Lorenz 
system given by eq. (3) with X, /*=0 (or eq.(5) with X=0), 
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direct benefit concerning the problem of "predicting predict¬ 
ability," which is currently a research topic of great interest. 

Finally, we remark that a similar chaos by chaos control¬ 
ling mechanism, having a different type of coupling, was 
recently proposed to be responsible for stabilization of the 
Earth’s obliquity by the Moon [Laskar et al., 1993]. 

Appendix 

Our controlling strategy results in the following dynamical 
system 


Figure 3. Lorenz attractor controlled by random trajectory of 
eq. (4): 1M1.3, 1.7), W*=0 and X=10.0. 


X =/(x) + k(y-x) 

y = g( y) + v(x-y)- 


while in Figure 3 the same attractor, modified by coupling 
the Lorenz system with eq. (4) with \=10.0, is presented. 
Although the attractor of Figure 3 is chaotic, its dimension 
is smaller than the dimension of the original chaotic attractor 
of Figure 2(a). Thus the dynamics of the modified system 
(5) is far more predictable than the dynamics of a single 
Lorenz model. Its Lyapunov dimension estimated from time 
series by Wolfs et al. [1985] methods is d h = 2.008. This 
dimension decrease is again produced by significant reduc¬ 
tion of positive Lyapunov exponent (X^0.69 in comparison 
with Xjssl.87 of the original attractor); 

We would like to note that, despite the fact that in neither 
of our numerical examples the conditions of Propositions 1 
and 2 have been precisely fulfilled, nevertheless application 
of our controlling procedure has allowed us to convert one 
type of chaotic behavior to chaotic behavior which is more 
predictable as we choose parameters X and /* close to those 
that fulfill the above mentioned conditions. For X and /x 
which are far away from the conditions of Propositions 1 
and 2 our controlling procedure does not work. 

In the cases considered in this letter the application of a 
continuous control scheme did not result in obtaining 
periodic behavior because the external periodic perturbation 
was not taken in the form of an unstable periodic orbit of 
original chaotic system. If that is done, a similar controlling 
scheme allows us to convert the original chaotic behavior 
into an appropriate periodic one [Pyragas, 1992, 1993, Qu 
et al., 1993]. It should be mentioned here that, if f=g and 
0 in eqs. (1), our controlling procedure simplifies the 
method of synchronization of chaos using continuous control 
[Kapitaniak, 1994]. 

In summary we propose that this procedure of controlling 
chaos by chaos can be treated as a possible mechanism for 
the so-called extended range atmospheric predictability 
observed in geophysical systems. The results of this simple 
coupling show the great potential influence of the behavior 
of chaotic system on that of another. It is also clear that, in 
a stationary state, when the variables x and y are close 
together, the control signals (2) are small (under conditions 
of Proposition 2 they converge to each other). This probably 
means that such a coupling might not allow for easy experi¬ 
mental verification in real geophysical systems. It is, 
however, a mechanism which may account for the existence 
of many unexpected short and relatively predicted features 
in the situations when strongly chaotic behavior might be 
expected. We hope that the results of this letter could be of 


where X, are real nonnegative parameters. Note that many 
systems can be put in the form of (Al) including the formula¬ 
tion [Smale, 1967] of the Turing reaction-diffusion theory 
[Turing, 1953] of morphogenesis or the evolution of two 
resistively coupled electrical circuits [Kapitaniak et al., 1993]. 

We assume that the dynamical systems 


y = g(y)> 

where x,yeR 3 and (A2) have chaotic attractors A,, A g and A 
respectively. Denote the projection of A on the subspace 
x =( x 1 ^ 2 , x 3 ) t by A x , and the subspace y=(y „y 2 j 3 ) T by A r 
Recalling the definition of topological equivalence of two 
chaotic attractors: namely an attractor A f is equivalent to 
attractor A g if there exists a homeomorphism h: R°-*R n such 
that h(A f )=A g we have the following propositions 

Proposition 1. 

(i) If f=g and 

1*0 = 0) -y(t = 0) I 

is sufficiently small, then there exists a value of &=X+/x, say 
k *, such that for k>k *, A x is equivalent to Ay. 

(ii) If f^g and X= oo , then A, is equivalent to Ag. 

(iii) If f¥=g and (x= then Ay is equivalent to Af. 

Proof: 

(i) First note that the inequalities 

I a if\ > L |*| 

4 = 1, i+j 

wherey=l,...,«, are sufficient for the stability of a matrix [a^ 
with negative diagonal elements. 

Denote w=x-y, so that from (Al) we have 

u = [-(* + |i )E*Ef\ vC ]u + 0{x,y) ^ 


= Au + O(x.y) 




1260 


BRINDLEY ET AL.: CONTROLLING CHAOS BY CHAOS IN GEOPHYSICAL SYSTEMS 


where Df is the Jacobian matrix of /, E is the unit matrix and 
0 (jc, y) represents the higher order terms. It is obvious that one 
can find k such that matrix is stable, that is u=0 is 

asymptotically stable, and x(?) approaches y(/) as t-*oo. Hence 
A x is equivalent to Ay (the homeomorphism h:R n ->R n is 
identity). 

(ii) Equation (Al) can be rewritten as 

ex = ef(x) + ( y-x ) 

y = g(y) + n(*-y), 

where e=l/\. If e=0, the last equation is equivalent to 
x = y 

y = gOO- 

Thus, A* is equivalent to Ag (again, the homeomorphism 
h:R n -*R n is identity). 

i 

(iii) The proof is similar as in the case (ii). 

The second and the third part of the Proposition 1 can be 
improved in the following way: 

Proposition 2. 

For sufficiently small 

I x(i = 0) + y(t = 0) I 

and e there exists t 0 such that x(t) converges uniformly to y(t) 
as e -* 0 + on all subsets of t 0 <t< oo. 

Proof: The proof is similar to the proof of Theorem 2 in 
[Kocarev and Kapitaniak, 1994]. 
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